The distance upon contact: Determination from roughness profile 
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The point at which two random rough surfaces make contact takes place at the contact of the 
highest asperities. The distance upon contact do in the limit of zero load has crucial importance 
for determination of dispersive forces. Using gold films as an example we demonstrate that for two 
parallel plates do is a function of the nominal size of the contact area L and give a simple expression 
for do(L) via the surface roughness characteristics. In the case of a sphere of fixed radius R and 
a plate the scale dependence manifests itself as an additional uncertainty Sd(L) in the separation, 
where the scale L is related with the separation d via the effective area of interaction L 2 ~ irRd. 
This uncertainty depends on the roughness of interacting bodies and disappears in the limit L — > oo. 

PACS numbers: 68.35.Ct, 68.35.Np, 12.20.Fv, 68.37.Ps, 85.85.+j 



I. INTRODUCTION 

The absolute distance separating two bodies is a pa- 
rameter of principal importance for the determination of 
dispersive forces (van the der Waals hj, Casimir [2|] or 
more general Casimir-Lifshitz force Q). The absolute 
distance becomes difficult to determine when the separa- 
tion gap approaches nanometer dimensions. This compli- 
cation originates from the presence of surface roughness, 
which manifests itself on the same scale. In fact, when 
the bodies are brought into gentle contact they are still 
separated by some distance do, which we call the distance 
upon contact due to surface roughness. 

We are interested in the dispersive forces when stronger 
chemical or capillary forces are eliminated. In this case 
do has a special significance for adhesion, which is mainly 
due to van der Waals forces across an extensive noncon- 
tact area [4]. The distance do is important for micro 
(nano) electro mechanical systems (MEMS) because stic- 
tion due to adhesion is the major failure mode in MEMS 
0| . Furthermore, the distance upon contact plays an im- 
portant role in contact mechanics p, is very significant 
for heat transfer [7[, contact resistivity [8|, lubrication, 
and sealing [9j]. In addition, it has also importance in the 
case of capillary forces and wetting [TO, TTTL fl2| , where 
knowledge of do provides further insight of how adsorbed 
water wets a rough surface. 

The distance upon contact do between a sphere and a 
plate pjj 03 plays a key role in modern precise measure- 
ments of the dispersion forces (see [HI for a review) where 
do is the main source of errors. In Casimir force measure- 
ments do is determined using electrostatic calibration. 
In this case the force dependence on the separation is 
known, and one can determine the absolute separation 
(see resent discussions [l6|, 03 Gil)- Even when the dis- 
tance is not counted from the point of contact [H, Qj], Qjl 
local realization of roughness as shown in this paper will 
contribute to uncertainty of the absolute separation. 

Independent attempts to define do were undertaken in 
experiments measuring the adhesion energy It was 



proposed [20[ to take do as the sum of the root mean 
square (rms) roughnesses of two surfaces upon contact. 
This definition is, however, restricted and can only be 
used for rough estimates as stressed in [20|. Obviously, 
the distance upon contact has to be defined by the high- 
est asperities. 

In this paper we propose a simple method for determi- 
nation of do from the roughness profiles of the two sur- 
faces coming into contact. For two plates it is explicitly 
demonstrated that do{L) is scale dependent, where L 2 
is the area of nominal contact. We discuss also applica- 
tion of our method to the sphere-plate configuration. In 
this case it is shown that do determined from the electro- 
static calibration can differ from that playing role in the 
dispersive force and the difference is scale (separation) 
dependent. 

In Sec. [Ill we report briefly the details of our film prepa- 
ration and characterization. In Sec. IIIII the roughness 
profiles in the plate-plate configuration are discussed and 
the main relation connecting do with the size of the nom- 
inal contact is deduced. The sphere-plate configuration 
is discussed in Sec. HVl together with uncertainty in do. 
Our conclusions are collected in Sec. [VI 



II. EXPERIMENTAL 

The surfaces we use in this study were gold films grown 
by thermal evaporation onto oxidized silicon wafers with 
thicknesses in the range 100 — 1600 nm and having dif- 
ferent rms roughnesses. A polysterene sphere (radius 
R = 50 /im), attached on a gold coated cantilever, was 
first plasma sputtered with gold for electrical contact, 
and then a 100 nm gold film grown on top of the initial 
coating. The deposited films were of uniform thickness 
and of isotropic surface morphology as was confirmed in- 
dependently with atomic force and scanning electron mi- 
croscopy on different locations. 

The surface profile was recorded with Veeco Multimode 
atomic force microscope (AFM) using Nanoscope V con- 
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FIG. 1: (Color online) AFM megascan of the 100 nm film (a) 
and the sphere (b). The insets show the highlighted areas at 
higher magnifications. 



troller. To analyze the effect of scale dependence, megas- 
cans of large area up to 40 x 40 /im 2 were made and 
recorded with the lateral resolution of 4096 x 4096 pix- 
els. The maximal area, which we have been able to scan 
on the sphere, was 8x8 am 2 (2048 x 2048 pixels). All 
images were flattened with linear filtering; for the sphere 
the parabolic filtering was used to exclude the effect of 
curvature. Figure [1] shows the images of the 100 nm film 
(a) and the sphere (b) on different scales. Approximately 
10 images of smaller size 500 x 500 nm 2 were recorded 
for each film and for the sphere to obtain the correlation 
length £ of the rough surfaces [21]. Finally, the electro- 
static calibration was used for the determination of the 
cantilever spring constant and do [22[. 



III. PLATE-PLATE CONTACT 

Consider first two parallel plates, which can come into 
contact. A plate surface can be described by a roughness 
profile hi(x,y) (i — 1,2 for body 1 or 2), where x and 
y are the lateral coordinates. The averaged value over 
large area of the profile is zero, (hi(x,y)) = 0. Then the 
local distance between the plates is 



d(x,y) =d - hi(x,y) - h 2 (x,y). 



(1) 



where d is the distance between the average planes. We 
can define the distance upon contact do as the largest 
distance d = do, for which d(x, y) becomes zero. 
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FIG. 2: (Color online) Distance upon contact as a function of 
the length scale. Dots with the error bars are the values calcu- 
lated from the megascans. The solid curve is the theoretical 
expectation according to Eq. J4). The inset demonstrates 
absence of the scale dependence for the rms roughness. 

It is well known from contact mechanics [23| that the 
contact of two elastic rough plates is equivalent to the 
contact of a rough hard plate and an elastic flat plate with 
an effective Young's modulus E and a Poisson ratio v. In 
this paper we analyze the contact in the limit of zero load 
when both bodies can be considered as hard. This limit is 
realized when only weak adhesion is possible, for which 
the dispersive forces are responsible. Strong adhesion 
due to chemical bonding or due to capillary forces is not 
considered here. This is not a principal restriction, but 
the case of strong adhesion has to be analyzed separately. 
Equation (pQ) shows that the profile of the effective rough 
body is given by 



h(x,y) = h±{x,y) + h 2 (x,y). 



(2) 



The latter means that h(x,y) is given by the combined 
image of the surfaces facing each other. 

Let Lo be the size of the combined image. Then, in 
order to obtain information on the scale L = Lo/2 n , we 
divide this image on 2 n subimages. For each subimage 
we find the highest point of the profile (local do), and 
average all these values. This procedure gives us do(L) 
and the corresponding statistical error. Megascans are 
very convenient for this purpose otherwise one has to 
collect many scans in different locations. 

For the 100 nm film above the 400 nm film the result of 
this procedure is shown in Fig. [2l We took the maximum 
area to be 10 x 10 fim 2 . The figure clearly demonstrates 
the dependence of do on the scale L although the errors 
appear to be significant. The inset shows the dependence 
of the rms roughness w on the length scale L. This de- 
pendence is absent in accordance with the expectations, 
while only the error bars increase when L is decreasing. 

To understand the dependence do(L) let us assume 
that the size L of the area of nominal contact is large 
in comparison with the correlation length, L ^> £. It 
means that this area can be divided into a large number 
N 2 = L 2 /^ 2 of cells. The height of each cell (asper- 
ity) can be considered as a random variable h [24]. The 
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probability to find h smaller than some value z can be 
presented in a general form 

P(z) = l-e"*W, (3) 

where the "phase" (j)(z) is a nonnegative and nondecreas- 
ing function of z. Note that ([3|) is just a convenient way 
to represent the data: instead of cumulative distributions 
P(z) we are using the phase (j){z). 

For a given asperity the probability to find its height 
above do is 1 — P(c?o), then within the area of nominal 
contact one asperity will be higher than do if 

e -*(*> (L 2 /0 = 1 or ^o)=ln(L 2 /C 2 )- ( 4 ) 

This condition can be considered as an equation for the 
asperity height because due to a sharp exponential be- 
havior the height is approximately equal to do- To solve 
(|4]) we have to know the function c/)(z), which can be 
found from the roughness profile. 

The cumulative distribution P(z) can be found from 
a roughness profile by counting pixels with the height 
below z. Then the "phase" can be calculated as <j)(z) = 
— ln(l — P). The results are presented in Fig. [3l It has to 
be noted that the function <j>(z) becomes more dispersive 
at large z. This effect was observed for all surfaces we 
investigated. To solve Eq. (jlj) we have to approximate 
the large z tail of <j)(z) by a smooth curve. Any way of 
the data smoothing is equally good, and our method is 
not relied on specific assumptions about the probability 
distribution. The procedure of solving Eq. is shown 
schematically in Fig. El and the solution itself is the red 
curve in Fig. [2l 

It has to be mentioned that the normal distribution 
fails to describe the data at large z. Other known distri- 
butions are not able satisfactory describe the data at all 
z. Asymptotically at large z the data can be reasonably 
well fit with the generalized extreme value distributions 
Gumbel or Weibull [25!, 26] . This fact becomes important 
if one has to know do for the size L, which is larger than 
the maximal scan size. In this case one has to extrapo- 
late cf)(z) to large z according to the chosen distribution. 
In this paper we are not doing extrapolation using only 
4>(z) extracted directly form the megascans. 

The observed dependence do(L) can be understood in- 
tuitively. The probability to have one high asperity is ex- 
ponentially small but the number of asperities increases 
with the area of nominal contact. Therefore, the larger 
the contact area, the higher probability to find a high 
feature within this area. 

Our result found in the limit of zero load will hold true 
if the elastic deformation of the highest asperity will be 
small (<C do). Applying Hertzian theory to an asperity 
of radius £/2 one finds the restriction on the load p: 

p « y^/9d~o (1 - v 2 ) _1 (dl/L 2 ) E. (5) 

If p = An/^Trdo is the van der Waals pressure (Ah is 
the Hamaker constant) then ([5]) for the Au parameters 
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FIG. 3: (Color online) Statistics of the surface roughness. 
Four 10 x 10 nm 2 images were used. The main graph shows 
the "phase" as a function of z. The red (light gray) curve 
is the best fit of the data at large z and the dashed lines 
demonstrate the solution of Eq. (j4}. The top inset presents 
the logarithm of the density function. The bottom inset shows 
the cumulative distribution. 

restricts do and L as (do/10 nm) 4,5 (L/10 /im) -2 ^> 0.3. 
This condition is true in the range of main interest. For 
the sphere-plane case (see below) Eq. (|5j) can be modified 
accordingly but in general the physical contact is not 
assumed for the sphere-plate configuration. 

IV. SPHERE-PLATE CONTACT 

The other question of great practical importance is the 
distance upon contact between a sphere and a plate. In 
the experiments [l3|, [3, EH, [22[ the sphere attached 
to a cantilever or an optical fibre approaches the plate. 
Assuming that the sphere is large, R^> d, the local dis- 
tance is 

d(x, y)=d+(x 2 + y 2 ) /2R - h(x, y), (6) 

where h(x,y) is the combined profile of the sphere and 
the plate. 

Again, do is the maximal <i, for which the local distance 
becomes zero. This definition gives 

do = max [h(x, y) - (x 2 + y 2 ) /2R] . (7) 

In contrast with the plate-plate configuration now do is 
a function of the sphere radius i?, but, of course, one can 
define the length scale Lr corresponding to this radius 
R (see below). 

As input data in Eq. Q we used the combined images 
of the sphere and different plates. The origin (x = 0, y = 
0) was chosen randomly in different positions and then do 
was calculated according to ((7|). We averaged do found 
in 80 different locations to get the values of c^ 771 , which 
are collected in Tab. H 

We can estimate the same value theoretically. A circle 
of a finite area L 2 is important in Eq. ([7|). Asperities 
of the size £ are distributed homogeneously within this 
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TABLE I: The parameters characterizing the sphere- film sys- 
tems (all in nm). The first five rows were determined from 
combined images (see text). The last row df} gives the values 
of do determined electrostatically. The last four rows were 
determined for R = 50 \im. 



circle. Then the averaged value of the second term in 
((7j) is L 2 /4ttR. The averaged maximal value of h(x,y) 
is the distance upon contact between two plates of the 
size L. This distance is the solution of Eq. (|4|). In this 
section we will denote it as cUq P (L) not to mix with do in 
the sphere-plate configuration. Then one can find do for 
the sphere-plate contact by maximizing Q on L: 

do = max [d p p (L) - L 2 /4irR] . (8) 

The solution of this equation defines d^ 1 and the scale 
Lr corresponding to the maximum. The values of d^ 1 
and Lr found from §8j) are given in Tab. [J for the radius 
R = 50 \im. 

One can see that d^ 1 is in agreement with d^ deter- 
mined from the combined images. Comparing it with 
the values d$ determined electrostatically one sees that 
in the first two columns the values of d$ are consider- 
ably larger. Moreover, the errors in d$ are smaller than 
in 4 m . 

We described do as the value determined from the area 
L 2 R and averaged over its different locations. Determina- 
tion of do from the electrostatic measurements did not 
undergo this type of averaging. As a result it is sensi- 
tive to the local roughness realization near the contact 
location. This explains why the errors in d$ are smaller: 
statistical variation of do from place to place is not in- 
cluded in the errors of d$ . 

Very different local values of do can be found and for 
this reason d$ can deviate significantly from the mean 
value. Choosing arbitrarily the contact locations in the 
image of the sphere and the 100 nm film we found, for ex- 
ample, that about 5% of the cases are in agreement with 
the measured value d$ = 17.7 ± 1.1 nm. One can imag- 
ine that the place of contact on the sphere has at least 
one asperity above the average. In the combined image 
the sphere dominates since it is rougher than the film, 
w sp h — 3.5 nm and w±oo = 1.5 nm. Because the sphere 
is rigidly fixed on the cantilever the same feature will be 
in the area of contact for any other location or other film. 
Already for the sphere above 400 nm film the high fea- 
ture on the sphere will not play significant role because 
the roughness of the film, 11*400 = 4.9 nm, is higher than 
that for the sphere. In this case we would expect that 
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FIG. 4: (Color online) Schematic explanation of additional 
uncertainty Sd in do (see text). The sphere in two positions 
is shown by the dashed (contact) and solid blue (dark gray) 
curves. The inset shows the variance of W as a function of 
the scale L or separation d. 



d$ has to be in agreement with the averaged value found 
from the image that is precisely what happens. 

Consider now the experimental situation when the dis- 
persive force is measured in the sphere-plate configura- 
tion. The system under consideration is equivalent to a 
smooth sphere above a combined rough profile h(x,y). 
The position of the average plane depends on the area 
of averaging L 2 especially for small scales L. The profile 
shown in Fig. 0] demonstrates different mean values in 
the left and right segments shown by the dashed black 
lines. Both of these values deviate from the middle line 
for the scale 2L (solid black line). The true average plane 
is defined for L — > 00. 

From Fig. [H one can see that do for L and 2L differ 
on Sd = do(L) — do(2L). To be more precise we can 
define the uncertainty in do as Sd(L) = do(L) — do, where 
we understand do as the value counted from the true 
average plane (L — > 00). The distance between bodies is 
then d = do + Sd(L) ± Ad, where Ad is the displacement 
from the contact point. The scale L is defined by the 
effective area of interaction L 2 = airRd (a = 2 for the 
electrostatic and a = 2/3 for the pure Casimir force). 
Suppose that do found from the electrostatic calibration 
can be considered as a true value (the electrostatic scale is 
large, L e \ — > 00) then in the dispersive force measurement 
the bodies are separated by d = do ± Sd(Ldi S ) + Ad with 
the related scale L^ is = \JcmRd. 

For a fixed L the uncertainty 5d is a random variable 
distributed roughly normally around Sd = 0. However, 
it has to be stressed that Sd manifests itself not as a 
statistical error but rather as a kind of a systematic error. 
This is because at a given lateral position of the sphere 
this uncertainty takes a fixed value. The variance of Sd 
is defined by the roughness statistics. It was calculated 
from the images and shown as inset in Fig. [U One has 
to remember that with a probability of 30% the value of 
Sd can be larger than that shown in Fig. [H 
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V. CONCLUSIONS 

In conclusion, it is shown that the distance upon con- 
tact depends on the lateral size of contacting plates and 
a simple formula describing do(L) is proposed. For the 
sphere and plate an additional uncertainty in the abso- 
lute separation d is revealed arising due to variation of 
the average plane position with the effective area of inter- 
action or equivalently with the separation. Its magnitude 
depends on the roughness of interacting bodies. 
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